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Q H ! Abstract 

The classical 0-C curves are discussed in different cases in which various period 
changes involved. Among them, the analytic 0-C curves with frequency, amplitude mod- 
ulations and with double modes are closely inspected, respectively. As a special, the 
light-time effect is illustrated. The features of period change noise and period change to 
metallicity are added at the end. 
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1 Introduction 



The 0-C diagram is a plot showing the observed times of maximum light (O) minus those 
calculated according to an adopted ephemeris(C) plotted as a function of time, mostly, the 
number of elapsed cycles. In the same way, the 0-C diagram can also be constructed by the 
difference between the observed times of maximum radial velocity and the times predicted from 
an adopted ephemeris. One may find the minima are used instead of maxima for some variables. 
In particular, the spectroscopic and photometric 0-C values are combined to produce a single 
0-C diagram. 

The employment of the 0-C diagram almost means normal or regular periodic light curve 
with a large amplitude is concerned. In other words, the times of maximum light can be 
determined sufficiently well from the observed light data. One can find '0' by fitting a single 
sinusoid with an assumed pulsation period to observations. '0' may be derived through local 
fit to the light curves around individual maxima as well. The latter is very useful in the case 
of asymmetrical light curves or light curves with changing amplitude such as the cycle-to-cycle 
or the day-by-day variations. 

The Calculated times of maximum light (C ra ) based on an adopted or estimated period 
(P es t) and an initial epoch or Time of maximum light (To). Assuming that the observed times 
of maximum light, (To, Ti, T2, • • • , Tv), are determined in different ways depending on different 
observers, we give 0-C series as 

(O - C) n = T n -(T + n P est ), n = 0, 1, 2, ,N (1) 

where n is the cycle number elapsed from the initial epoch at the observed point. Followings 
refer to different cases of period change in the point of O-C curve. 

*published in 1999, Publ. Beijing Astron. Obs., 33, p.17 
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2 Different O-C Curves 

2.1 Period Changing Nothing 

At the beginning, we consider the case of period P keeps constant over the time span of 
observations covered. By assuming that 

P est = P + 5P, T cst = T + 5T (2) 

then we can compute O-C as 

(O - C) n = T n - [T + 5T + n(P + SP)} (3) 

If 5P=0, namely 

(O - C) n = T n - (T + 5T + nP), ((O - C)) = -5T (4) 

the O-C diagram will consist of points scattered around a horizontal line: 

(O - C) n = -5T (5) 

by assuming T n — (T + nP) = O(0. 0), the same below. 
If 5T =0, that is 

(O-C) n = T n -(T + nP)-n5P (6) 
Clearly {O — C) n decreases with n increases, the O-C plot will be a straight line: 

(0-C) n = -n5P (7) 

If 5P > 0, that is P est is longer than P, the O-C line will slop downwards; otherwise, 5P < 0, 
slop upwards. A vertical shift of the (O — C) n line relative to 5P=0, 5T =0 lines in Fig. 1 tells 
us a shift in the epoch, the change in the phase of variation in other words. 

2.2 Period Changing with Time 

Now, if the period is changing with time, may be increasing or decreasing. Let P(t) stands for 
a changing period with a slow, constant rate (3 as caused by stellar evolution, we write 

P(t) = P + (3t (8) 

The derivative of the equation above gives the rate of period changing = (3. The stellar 
radius or mass changes due to evolutionary effects or mass loss often results in a slowly changing 
in period. We have 

C n = T + nP , O n = T + nP(T n ) , . 

(O - C) n = n[P(T n ) - P ] = n(3T n , n = 0, 1, 2, ,N W 

By denoting AP as the mean change during one period, when T x = P,T 2 = 2P, , T N = 

NP, we obtain 

P(T 1 )=P + AP, , P(Tv) = P + iVAP (10) 
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Figure 1: The 0-C diagram as a constant period. The slop line, (O — C) n = n5P, refer to the 
estimated period P cst = P + 5P assuming the intrinsic period P is constant. Here, ( s >P=0(P)<^ 
P. Dashed lines, (O — C) n = ST , correspond to the cases in which a correct period value 
is employed but the initial epoch has an uncertainty of 5T . Note that the mirror parts are 
omitted. 
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Figure 2: The O-C diagram as the period is changing slowly with time. (O — C) n = APn 2 . 
The solid curve for AP > 0, the dashed curve for AP < 0. 
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So during the time span of the observations concerned, the change of the period is 

(O - C) n = n[P(T n ) - P ] = APn 2 (11) 

This is a parabola equation (see Fig. 2). 
With assumption that 

C no = T + P n , O no = T (n ) + P n (12) 

we have 

(O - C) n = (O - C) n „ no + (0- C) no = AP(n - n ) 2 + (T (n ) - T ) (13) 

Now the equation demonstrates the information of both the epoch T and the period change 
AP are included in the O-C diagram. The measured accuracy of period changes increases as 
n 2 . Hence it is more important to extend the observations for a well-studied variable star. 

The predicted evolutionary period changes of variable stars is commonly difficult to be 
revealed in the O-C diagram because the evolutionary changes may still be too small to be 
detected after a few tens of years. In that cases, the O-C diagram should be essentially a 
straight line. In special, in those stars evolving to the red toward the end of their horizontal 
branch lifetimes, the evolutionary period changes, while still small, ought to be large enough 
to be detected. In these cases, the O-C diagram should be well represented by a parabolic 
curve, which should indicate a slow rate of period increase. A few stars might be caught in 
the interval of core instability, and could show period increases or decreases much larger than 
dP/dt— 0.025 days per million years. 

Prager (1939) noted that three field RR Lyrae stars' periods did not change at a constant 
rate, but were better modeled assuming abrupt period changes. Abruptness, of course, does 
not accord with a model in which period change results from smooth evolution. Usually, the 
period changes revealed in the O-C diagram were more abrupt than a perfect parabola would 
allow. The period changes are small compared to the length of the period itself but because 
many cycles have passed by, the cumulative effects on the O-C diagram may be large. In some 
stars, the observations can be well fit by assuming that the period is changing at a constant rate 
(parabola O-C diagram), but many others follow Prager's stars seem to have abrupt changes 
in period. When fit by a parabola, the derived rates of period change are usually in the range 
of —1 to +1 days per million years, or a change of less than 10 seconds during 50-100 years in 
which the stars have been under observation. 

2.3 The Phase-shift Diagram 

In addition to the O-C diagram, an analogy, the phase-shift diagram is also applied to interpret 
period changes in light (e.g. Percy et al. 1980, Breger 1990). Followings show the phase-shift 
diagram functions exactly the same as the classical O-C diagram. We define quantity /, as 
below: 

/. - < 14 > 

-Test 

here ti is the observed times such as times of maximum light of pulsating stars, to is any 
arbitrary reference time, P est is an estimated period. The fraction remained of fi constructs 
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the phase fa, and the integer of /j produces the cycle number N. A plot of the phase against 
the cycle number is thus called the phase-shift diagram. Under the assumption that the period 
P is changing uniformly at a rate of (3 described by eq.(8) where Po is the period value at t , t 
is the time interval passed by from to, hence 

P i = P + P(t i -t ), i = 1,2,3, ,N (15) 

If we take the mean of Po and Pj, the period at middle point of the two periods is 

(P) = \(Pi + P o) =Po + \m - to) (16) 

we have cycle numbers from P to Pf 

U — to ti — to n Po 

71 = /n\ = n > \ot+ 7T> (*» ~ *0 



(P) P + l(3( ti -toy Vl u; l-i/3n 
Such that we find fa: 

, U — to nP , 

* , = -RT-" = p-.(i-i/J n )- B (17) 

Without loss of generality < 1, we make use of series expansion that 
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T nT- = E(oH n «i + > (is) 

1 n=0 Z Z 

therefore, 

0„ = ^(1 + Un) - n = (Po - P est )^- + = (P - Pest)^ + ^/3n 2 (19) 

-Test 1 J^est ^-r est -Test ^ 

(P /P cs t ~ 1, n — 1, 2, 3, , N). The equation presents a parabola about cycle number n. 

Quantity /3, is generally used to characterize the period change rates of variables. It is the rate 
of period change determined through fitting a parabola to the O-C diagram, and is normally 
presented in the unit of days/days or days/year. A unit of days per million years for /3 is chosen 
by some authors. -p^£ is sometimes presented as the period change rate in unit cycle/year. 



3 Analytic O-C Curves 

We describe the light brightness of a variable star with a sinusoidal waveform as: 

L(t) = Lo + A sin(27rf/P) (20) 

where P is the period of light variation, A is the amplitude of the sinusoidal wave, L is the 
zero point. It is well known that the times of maximum light occur at L'(t) = 0, which 

yields O n = {T ,T 1 ,T 2 , ,P/v}- At the same time, sin — > 1 and cos 2lT p n — ► 0, at 

the times of maximum light. The period changes in the light curves are discussed in several 
different sources below. 
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3.1 Frequency Modulation 

For frequency or period modulation, theoretical light curve is defined by 

L t (t) = L Q + A sm(2irt/P(t)) (21) 

L' { (t)=0 produces O n as 

Pit) 

O n = ^(2n + l), n = 0,l,2, ,N (22) 

C n are naturally derived from L(t) = L + A sm(2irt/P), so 

(O - C)„ = \(2n+ l)(P(t) -P), n = 0, 1, 2, ,N (23) 

Recalling eq.(8): P(t) = P + fit = P + j3Pn which means (5 is the period change over one cycle, 
thus 

(O - C) n = \fiPn 2 + \fiPn (24) 
this is a parabola with respect to cycle number n. 

3.2 Amplitude Modulation 

The change in period caused by the variability of amplitude called amplitude modulation is 
described as: 

Ljf) = L + A(t) sin(27rt/P), A(t) = a sin(27rf/P A ) (25) 

A(t) refers to a sinusoidal amplitude modulation which causes a secondary period Pa in the 
light curve. The calculus of L a (t), L' a (t), gives observed times of maximum light O n : 

1 27rO„ . 2ixO n 1 . 2vrO„ 2vrO n 

— cos sin sin cos = 26 

Pa Pa P P Pa P K J 

By making use of the replacement of O n = (0 — C) n + C n just for the triangle functions involving 
P and expanding them. Without loss of generality we presume 

(0 _ C) „ <<P; ^M^k.M^k (27) 

and with the help of the features mentioned in the beginning of this section: sin 2lT p" — >■ 1, 
we find 

1 2ixO n .2n(0 - C) n 2ixC n 2ix{0-C) n , 
- cos — ( cos — + cos ) + 

1 . 2nO 2tt(Q - C) n 2nC n 2n(Q - C% 

P Sm ~^ (C ° S P C ° S ~P P } = ( } 

Expanding similarly now the two functions including Pa, and recalling that cos ^pp- — > 0, 
cos M °~ c)n — > 1, L' a (t) = becomes 

1 27r(0-C)„ 27rC n 1 2Tx(0-C) n 
cos cos ■ 



-Fa -Pa -Pa Pa Pa 



1 [2n(0 - C) n ] 2 2irC n 1 2n(0 - C) n 2tt(0 - C) n n 

— — — cos — — cos = 

P PP A P a P P P a 

Note that the coefficient of the third term of the equation above is infinitely small at higher 
order of ^°~^ n which is set to almost zero as our original consideration. This yields 

2tt(0 - C)n(~p2 - p2 COS Y K ^ ~P~ K Ta ~Ta~ = ( ^ 

When Pa is so large enough that < 1, we can make following simplifications: 

2x(Q-0. h1 _MQ-0. [2«{Q-C) n F 

Pa Pa Pa K ' 



Consequently, we get 



p A cos 2w< '" 



i0 ~ C)n= 2^TT{W^W (31) 

This relation shows the shape of 0-C is predominated by the modulation period Pa- If Pa ~ P, 
like cycle-to-cycle modulation, 

that indicates a shift of Pa/2ti with respect to the zero line showing constant periodic light. 
In the case of great difference between Pa and P, if Pa 3> P, the long-term amplitude modu- 
lation, stimulation shows only a tiny scatter around the zero of the O-C, while Pa <C P, the 
high frequency amplitude modulation, the feature of the O-C is of the behaviour of periodic 
amplitude variation, we have after neglecting the term including P in eq.(29), 

fn ^ 2tt(0 - C) n 2nC n 

(O -C) n = — cos cos — — (32) 

An F A fA 

A sinusoidal light curve with an amplitude modulation may have a significant variation in O-C, 
while a saw-tooth light curve with the same amplitude modulation will show an absolutely flat 
O-C. 



3.3 Double Mode Case 

As for double-mode variables with small amplitudes, we choose: 

L d (t) = L + Asm-— + Bsin{-— + 0) 
Pi P2 



Let L' d (t)—0 for finding O n and so (0 — C) n . Replacing O n with (O n — C n )+C n = (0 — C) n + C n . 

Q-C)n I • 2n(Q~C) n _ 2tt(Q-( 



With assumption of (O - C) n < P u so cos 2 ^° p c ^ — > 1, sin 2 ^°~ c ^ « *<o-c) n and 



sin 2n (° c *) n ~ 2 Ao c) n ^ and reca }} that we wan t to find P 2 in the O-C diagram. Therefore, 
the expanding of L' d (t)=0 is equivalent to 

P 2 A 2n(0-C) n . 2nC n 2vr(0 - C) n . .2nC n P 2 A 2irC n ,2nC n 

— — — sm — sm — h<z> = — • — -cos — hcos — h<z> 

Pi B P l P 1 P 2 v P 2 ; Pi B P 1 v P 2 YJ 
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thus we have (O — C) n like 

1 ^|cos^ + cos(^ + 0) 
( ° ~ C)n= 2i Jp sin ^ + J sin(% + 0) (33) 

At the times of maximum light, we know cos — > 0, sin — > 1. Finally, we get 

P 2 cos(^p + 0) 

( °- c) " = Mapf + si ;(^ + « (34) 

For the two different periods, P2/P1 ~ 0.78 is the normal case for the large-amplitude S Set 
stars. Taking a common amplitude ratio value of 0.30 and a second period of 0.03 d we may 
simplify eq.(34) to be 



(O - C) n 



0.005 cos(^ + 0) 
0.2 + sin(^ + (f) 



Beating phenomenon can also be involved in the feature of the O-C curve when the two 
periods are closely spaced with near amplitudes. Furthermore, the effects of the beating will 
be particularly obvious in the times of maximum light, giving a large O-C amplitude. An 
interference of two oscillations with the same amplitude A for a suspect (5 Cep star, 53 Psc 
(Wolf 1987), leads to the amplitude modulation of the form: 



. / / 111 

A h = A^(l + cos-), ^ = ^"^ (35) 

here Pb is the beat period. Based on this equation, the time required for the amplitude to reach 
a certain value A Q b s , counted from the epoch of Ab=0, can be derived as 

At » ^ (36) 

One more example is the light-time effect in a binary system if P 2 ^$> P±, from which the O-C 
curve may be interpreted, 

^"-sk-!" 8 ^*) (37) 

The information revealed by the O-C diagram depends on both 'O' constructed from max- 
ima (often used for pulsating stars) or minima (as in binary) and the methods used to derive 
the 'O' from the data. Usually four ways to find 'O', we take 

• the time of maximum brightness for each single observation; 

• the time of maximum brightness in a smoothed light curve; 

• the phase of the maximum of a representative mean light curve fitted to all or a portion 
of the observations over that cycle; 

• the time of maximum derived from an exact solution for a Fourier series which has been 
fitted to the data. 
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4 The Light-Time Effect of A Binary 



If a pulsating star is a member of a binary system, an independent estimate of the mass from 
stellar pulsation theory can be given. This is the intention of the detection of the light-time 
effect (LTE) of a pulsating star. 

4.1 Examples of Three (3 Cep Stars 

Pigulski has throughly investigated the nature of the period changes of three large-amplitude (3 
Cephei-type stars: f3 Cephei itself (the prototype of this class), a Scorpii and BW Vulpeculae 
in terms of the combined 0-C diagram constructed from both the light and radial velocity 
data. It was shown that the observed period changes were caused by the superposition of a 
constant-rate period change originated from normal stellar evolution and of the light-time effect 
in a binary system. 

Pigulski and Boratyn (1992) proposed that the observed changes of (3 Cephei (=HR 8238=HD 
205021, V=3.23, Bl IV) can be understood in terms of orbital motion of the variable in a binary 
system, in which the secondary is the less massive companion resolved by speckle interferom- 
etry. In fact, the companion in the system of (3 Cep was found in the speckle interferometric 
way by Gezari et al. (1972). By assuming that the orbital motion is the unique cause of the 
observed variation of the pulsation period, 

Ppuis = P (l + ^f) (38) 

Kad is fully described by six elements of the spectroscopic orbit: 7, radial velocity of the 
system's barycentre; P or b, orbital period; Ki, semi-amplitude of the radial velocity curve, due 
to orbital motion of the variable; e, eccentricity of the orbit; u>, longitude of periastron; T , 
the time of periastron passage. The analyses of a single O-C diagram constructed in terms of 
both the photometric and radial velocity times of maximum offered by Pigulski and Boratyn 
(1992) show that the apparent changes of the pulsation period of f3 Cep can be explained by 
the light-time effect, induced by the orbital motion of the variable around the barycentre of the 
variable-speckle companion system. 

a Scorpii (=HR 6084=HD 147165, Bl III) is a quadruple system. Its light variability can 
be described as a superposition of four periodic terms( Jerzykiewicz & Sterken 1984), of which 
the main pulsation period P 2 = d . 24684 and Pi=0 d . 23967 are also observed in radial velocity. 
The combined O-C diagram using the times of maximum radial velocity and light displays 
both an increase and a decrease of the main pulsation period. The period changes seen in the 
O-C curves can be explained neither by an evolutionary effect nor by LTE alone. 

(O - C)lt = (0-C) + d .000015 E - l d .29 x 10~ 10 E 2 (39) 

In fact, both effects contributed to the observed changes of the main pulsation period of a Sco. 
And LTE is caused by the speckle tertiary. Due to the observations do not cover the whole 
orbital period of the tertiary, the (O — C)lt diagram cannot help to give a more accurate value 
of the orbital period than that provided by Evans et al. (1986), is of the order of 100-350 yr. 
In this sense, the spectroscopic elements are unavailable from the (O — C)lt diagram. But it 
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is sure that the changes seen in the (O — C)lt diagram can be fully explained by the orbital 
motion (Pigulski 1992). 

BW Vulpeculae (=HR 8007=HD , V=6.55, B2 III) is a monoperiodic large-amplitude (3 
Cep star. Odell (1984) and Jiang Shi- Yang (1985) showed that the parabolic fit to the O-C 
diagram of BW Vul yields periodic residuals. Both authors attributed these residuals to LTE 
induced by the motion of the star in a binary system with an orbital period of about 25 yr. 
Pigulski (1993) showed that the combined effect of the increase of period with a rate of +2.34 
s/cen and LTE in a binary system (P orb =33.5±0.4 yr) fully accounts for the observed changes 
of the pulsation period of BW Vul. The light-time effect in the (O — C)lt diagram is, for the 
time being, the only observational evidence for the presence of a companion to BW Vul (see 
Fig. 3). For BW Vul, the times of minimum light and of minimum radial velocity were used in 
the determination of period changes. One ephemeris of minimum radial velocity was adopted 
so as to produce a single O-C diagram. 

All periodic variables with stable periods such as pulsating stars, eclipsing binaries and 
pulsars offer the possibility of detecting unseen companion(s) by means of the light-time effect. 
A pulsating star is a promising candidate for the detection of LTE provided that it belongs to 
a binary or multiple system. A simplest case is the LTE in a double system with a circular 
orbit. In the O-C diagram, the changes of the period P of a signal emitted by the visible 
component (pulsating primary) due to LTE will have the form of a sine curve with semi- 
amplitude equal to P or b K 1 /2%c, where 2K 1 is the range of the primary's radial velocity curve, 
P or b is the orbital period. The amplitude of LTE depends proportionally on P or b and K\ and is 
equal to the time which light needs for passing the projected orbit. For a given orbital period, 
short-period variables are more suitable for the study of LTE than the long-period ones because 
the visibility of LTE depends on the P or b/P ratio. Pigulski (1995) listed 10 objects including 
eclipsing binaries for which period changes attributable to LTE are observed. 

4.2 The Binary Model of The 5 Set Star CY Aqr 

The possibility that CY Aqr is a spectroscopic binary (see Fig. 4) immediately leads to the hope 
that a mass can be determined. In Figure 5 we show the relationship between the masses of the 
two companions for assumed values of the orbital inclination. The curve for % = 90° sets a lower 
limit to m 2 at each m\. Recent results given by McNamara et al.(1996) indicate: mi=l.O6M , 
M v —2.5. For this mass value of the primary star m 2 rounds 0.15 M Q . The companion must 
be fainter than CY Aqr by about 8.8 mag. We stress that, the binary model seems to us to 
be the most plausible interpretation of the variable times of maximum light. It remains only 
a reasonable hypothesis until new confirming data are available. The simplest and most easily 
accomplished test is to examine the radial velocity of CY Aqr with a 1.4 km s _1 amplitude 
predicted by the model. 

5 Period Change Noise 

Some stars appear to be unchanged in period over the time span of observations; others have 
undergone multiple changes in period. Some stars have exhibited both an increase and a 
decrease in period, as can be seen in the O-C diagram of RR Gem in Fig. 6. Stellar evolution 
is not expected as we know to act so quickly as to cause period changes in both directions in 
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Figure 3: The 0-C diagrams showing the light-time effect. Upper panel: (O — C)ur for the 
j3 Cep star BW Vul, the changes seen in this figure are interpreted as the light-time effect 
only. Continuous line represents the best fit, from which spectroscopic elements of the orbit are 
determined. After Pigulski (1993). Lower panels: sinusoidal variations in the O-C curves of the 
Cepheids FN Aql and RX Aur. Open circles denote the residuals obtained from photographic 
observations; photoelectric measurements are marked with dots. The size of these symbols 
refers to the weight assigned. Courtesy of Szabados (1988). 
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Figure 4: The 0-C diagram of CY Aqr (cycles) and the curve fitted using both a parabolic 
curve and a trigonometric function (solid line). 
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Figure 5: Solutions to the mass function for various choices of sini . 
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Figure 6: The O-C diagram for the times of maximum light of RR Gem, after Smith (1995). 

so short a span. Moreover, the sizes of the observed period changes are often much larger than 
expected from theory. We call this difference as "period change noise". In this regard, it is 
not unquestionable to attribute the period changes observed in any individual star to evolution 
alone. 

The period change noise is a result of random mixing events associated with the semiconvec- 
tive zone of the stellar core (Sweigart and Renzini 1979). The changes in the internal structure 
of a variable star arise in two ways: through the gradual change in core composition caused by 
nuclear burning, but also through a composition redistribution in the deep interior caused by 
convective overshooting and the formation of a semiconvective zone. Whereas this composition 
redistribution is assumed in standard calculations to proceed smoothly on the same timescales 
as nuclear burning. The nuclear reaction and convection need not be perfectly coupled at every 
instant of time. For an actual star, the composition redistribution is a discontinuous process. 
There might then be many mixing events, each slightly altering the internal of the star, and 
each leading to a change in period. Furthermore, different types of mixing events can account 
for both period increasing and decreasing. The observed timescales of period changes are in at 
least rough agreement with the mixing event hypothesis. Nuclear burning would be expected to 
produce a slow rate of period change which would be essentially unvarying over the span of a 
century. The mixing events, on the other hand, if they took place quickly would produce period 
changes which would seem abrupt; to the contrary, on a longer timescale they might produce 
period changes which, for one hundred years, would seem constant in rate. Of course, over a 
sufficiently long period of time, the period changes due to nuclear burning and discrete mixing 
events at the semiconvective zone must average to the value predicted by those evolutionary 
calculations which assume that the composition redistribution takes place smoothly. 

As another origin for the period change noise, proposed by Stothers (1980), which invoking 
hydromagnetic effects (such as for RR Lyrae stars). It is a common phenomenon involving 
stellar magnetism. 
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Third, Dearborn et al. (1990) suggested that exotic particles, such as weakly interacting 
massive particles (WIMPS), might produce changes in the pulsation periods (say of RR Lyrae 
stars). The existence of these particles has been postulated to explain the apparent existence of 
large amounts of dark matter in galaxies. These particles might become trapped within stars, 
where they can provide a new mechanism of energy transfer. In the case of horizontal branch 
stars, WIMPS might cause thermal pulses to occur on the Kelvin-Helmholtz timescale for the 
stellar core. These thermal pulses would produce period changes. The timescale calculated for 
these changes, however, on the order of 10 6 yr, is too long to explain the shorter term period 
noise. 

Fourth, the light travel time effect because of orbital motion in a binary stars system may 
be another possible source for apparent period changes as in RR Lyrae variables identified by 
Coutts (1971), and in f3 Cephei stars inspected by Pigulski & Boratyn (1992) and Pigulski (1992, 
1993). The motion of a star in a binary system could introduce a long term periodic oscillation 
in the O-C diagram for the time of maximum light, as the star alternately approaches and 
recedes from us. It is difficult, however, to build a case for the binary nature of any variables 
from O-C data alone. The examples are comparatively few. Saha and White (1990) have 
argued on the basis of O-C data and radial velocity observations that TU UMa, a RR Lyrae 
stars, may be a binary; Barnes III and Moffett(1975) and Moffett et al. (1987, 1988) evidenced 
the duplicity of SZ Lyncis, a dwarf Cepheid; others as BS Aqr (Fu Jiang-Ning et al. 1997), CY 
Aqr(Zhou Ai-Ying & Fu Jiang-Ning 1998), BE Lyn (Kiss & Szatmary 1995, Liu et al. 1994, 
Tang et al. 1992, Liu et al. 1991), BW Vul (Pigulski 1993) and so on (we refer the reader 
to examples preceded by Pigulski 1993, 1995). However, does the binary hypothesis by itself 
provide an adequate explanation for the period change behavior? 

Finally, mass loss might alter the theoretical period changes. For stars on the horizontal 
branch lost their masses when they were in the instability strip, that is, when they were RR 
Lyrae stars. But, the predicted rates of period changes were not significantly changed for mass 
loss rates as high as 10~ 9 M o /yr. 



6 Period Change to Metallicity 

As a variable star evolves through the H-R diagram, crossing its various stages of evolution, the 
metal abundance or chemical composition changes as well as its period does. In other words, 
at each given effective temperature periods vary with metallicity. Thus it is natural there is 
some correlation between the two observable quantities. In real, the period change depends on 
the change of chemical abundance. We write the relationship of the period change with respect 
to metallicity, P — [Fe/H], as 

The spectral type determined from the hydrogen Balmer lines (termed Sp(H)) was often 
later than the spectral type determined from the prominent Ca II K line at 3933 A (termed 
Sp(K)), which investigated by Preston (1959) under a low resolution spectroscopic survey of 
RR Lyrae stars in the solar vicinity. On the basis of the difference, Preston then proposed AS 
index as an indicator of metal abundance (hereafter the AS method): 

AS = 10[Sp(H) - Sp{K)\ (40) 
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which was later on confirmed by the analyses of high resolution spectra of RR Lyrae stars 
(Preston 1961; Butler 1975) and of model stellar spectra (Manduca 1981). Butler(1975) derived 

[Fe/H] = -0.23 - 0.16AS (41) 

Blanco (1992) further improved above as 

[Fe/H] = -(0.34 ± 0.02) - (0.18 ± 0.05)AS (42) 

Based on the globular cluster metallicity scale (adopted by Zinn & West 1984, and Zinn 1985), 
Suntzeff et al. (1991) got following calibration, 

[Fe/H] = -0.408 - 0.158 AS (43) 

is about 0.2 dex more metal-poor than Butler's scale. In the case of RRab stars, lower AS 
(high metallicity), smaller period changes; higher AS (low metallicity), bigger period changes. 
Iron abundance [Fe/H] is inferred from the strength of the calcium H-line. Certainly, such 
a calibration does not require that [Fe/H] = [Ca/H]. Thus, it does require that there exist a 
one-to-one correlation between iron and calcium abundances. In fact, [Ca/H] appears to be 
slightly high compared to [Fe/H] for metal-poor stars in general. Manduca (1981) provided an 
approximate relation as 

[Ca/H] = 0.B[Fe/H] (44) 

for theoretical calibration of the AS index. 

In addition to the method, other metallicity indicators have also been developed. The 
photometric indices among which are Sturch's(1966) broadband ultraviolet blanketing index 
5(U-B), the photometric k-line index (k — b) 2 (Jones 1971, 1973), the Stromgren mi index 
(Epstein 1969, Epstein and Epstein 1973), and the Walraven system 5[B-L] index (Lub 1977, 
1979). In brief, metallicity can be estimated from A«S(Preston), (k — 6)2(Jones), A(B — L) 
(Walraven), 5(U — B) (Sturch) and 8mi (Stromgren photometry). For the range 0< A S<12, 
we have relationships below among the approaches: 

AS = -845(U - B) + 14 

AS = -44(fc - b) 2 + 10 (45) 
AS= -99(toi) + 16 

(Butler 1975) and 

AS = 11.9 - 194A[£ — L] + 815(A[£ - L]f (46) 

(Blanco 1992). If metallicity did not change over the pulsation cycle of the star, these indices 
would remain constant. Rodriguez et al.(1990) failed a test of this internal consistency in the 
fact of SX Phe and the large amplitude 5 Set stars. Instead, they concluded that the m\ index 
variation is larger when the metal abundance is smaller. 

7 Ending Remarks 

The O-C diagram provides us the accurate period and maximum epoch of a variable star, a 
detector of its period and phase changes, and the information about the spectroscopic orbital 
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elements if it is a binary. Plus, the 0-C data provide constraints on stellar evolution theory. 
For binary systems, O-C analyses of eclipses can reveal very gradual changes due to tidal 
interactions, or even to gravitational radiation. 

The information of period changes is also presented as a phase diagram instead of a classical 
O-C diagram. Various shapes of the O-C curves correspond to following different cases: 

1. The O-C diagram for a star without measurable changes in period is a straight line. 

2. If the period of the star is constant, and if the correct period has been adopted, points 
on the O-C diagram will scatter around a straight horizontal line (see Fig. 1). 

3. If the period of the star is constant, but the adopted period is too long or too short, the 
straight line will slop up or down (see Fig. 1). 

4. If the period of the star is changing at a slow, constant rate (i.e. P(t) = P + j3t where (5 is 
small), then a good approximation of the O-C diagram can be represented by a parabola 
(see Fig. 2). 

5. If the period change is caused by the light-time effect in a binary system, the light curve 
appears to be sinusoidal (see Fig. 3). 

6. Broken lines — abrupt or jump changes of period (see Fig. 6). 

In order to interpret an O-C diagram correctly, it is necessary to know a number of cycles 
have elapsed between two observed maxima or since an initial epoch. The determination of 
elapsed cycles from an epoch is not always easy and reliable especially while the period of the 
star has been changing significantly. 

All common period-finding techniques assume that the periods are constant over the time 
span of the observations and that the phases are stable. However, these conditions are not 
always met by variables. Therefore, the O-C method is, in fact, one of the necessities of period 
analysis. A short review on the O-C method we refer to Willson (1986), a detailed theoretical 
study may be also referred to the paper series published in MNRAS (Koen &Lombard 1993, 
1995; Lombard & Koen 1993; Koen 1996). Koen and Lombard had paid scrupulous attention 
to the O-C methodology. They pointed out that the problem suffered from within the O-C 
method that the series of intervals among times is non-stationary so that the applications of 
this technique may lead to spurious conclusions. As well, they modelled the period changes 
in variable stars on a quantitative basis, especially supplied tests for the presence of intrinsic 
scatter in the periods of the long-period pulsating stars. 
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